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Abstract. Mass determinations are difficult to obtain and still fre- 
z^» \ quently characterised by deceptively large uncertainties. We review below 

the various mass estimators used for star clusters of all ages and lumi- 
^ \ nosities. We highlight a few recent results related to (i) very massive old 

star clusters, (ii) the differences and similarities between star clusters and 
CO \ cores of dwarf elliptical galaxies, and (iii) the possible strong biases on 

mass determination induced by tidal effects. 

> 

CO 

o 
o 



23 



Ctf 



X 



1. Introduction 



Open and globular clusters are respectively located, because of their local def- 
inition, in the plane and in the halo of our Galaxy. This review will not make 
Qh' any formal distinction between these two kinds of star clusters since (i) the 

above definitions apply only to our Galaxy and (ii) we do not know if there 
are any genuine differences in their formation mechanisms. Consequently, we 
may use the words "open" and "globular", but we shall essentially mean "star 
clusters", whether they are light or massive, or young or old. This approach 
is also justified by the fact that, in simulations of the Galactic globular cluster 
system, the dynamical evolution of an initial mass distribution, of either Gaus- 
sian or power-law type, leads always to a predicted distribution consistent with 
observations: light star clusters dissolve rather quickly, while heavy ones survive 
longer (Baumgardt 1998, 2001). See also Zhang & Fall (1999) in the case of the 
populations of star clusters of the Antennae Galaxies. 

Individual masses of star clusters are not easy to measure. As an example, 
we shall quote the four different mass estimates given for the giant Galactic 
globular cluster wCentauri by Ogorodnikov etal. (1976). From a few stellar 
radial velocities Mtot — 7 x lO 5 A / f0, from the same radial velocities corrected 
for the effect of the internal rotation of the cluster M.tot — 3 x 10 6 A4q, from low- 
quality proper motions measured on photographic plates M.tot — 3 x 10 7 \M®, 
and from gravitational focusing from the same proper motions Mtot — 2.8 x 
10 8 M.Q. These four different approaches provide results differing by 3 orders of 
magnitude! Fortunately, things have improved since then, but even the best mass 
determinations of star clusters remain rather uncertain, typically by a factor of 2. 
It is interesting to realize that the age of the Universe, intuitively more difficult 
to determine than the mass of a nearby star cluster, is more accurately known 
since the Hubble constant H is now estimated at a better than 20 % level. 
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2. The Giant Galactic Clusters lo Centauri and 47Tucanae 

Significant improvements in the quality and numbers of radial velocities have 
allowed the determinations of reliable mass estimates through the use of various 
dynamical models. On the one hand, there is the parametric approach. For 
example, Gunn & Griffin (1979) developed multi-mass models whose distribution 
functions / depend on the stellar energy per unit mass e and the specific angular 
momentum /. Such models are spherical and have a radial anisotropic velocity 
dispersion (v% ^ v | = vV). They are called King-Michie models and associate 
the lowered Maxwellian of the King model with the anisotropy factor of the 
Eddington models: 

f(e, I) oc [exp(-2j 2 e) - exp(-2j 2 £t )] exp(-j 2 / 2 /r 2 ) (1) 

On the other hand, one never knows which of the assumptions underlying such 
a model are adhered to by the real system and which are not (Dejonghe <fe 
Merritt 1992). These arguments suggest that it might be profitable to interpret 
kinematical data from globular clusters in an entirely different manner, placing 
much stronger demands on the data and making fewer ad hoc assumptions about 
the distribution function / as well as the gravitational potential <&. Ideally, the 
unknown functions should be generated non-parametrically from the data, in an 
approach pioneered by Merritt (1993, 1996). We provide hereafter the results of 
two studies, parametric and non-parametric, respectively, of the globular cluster 
(j Centauri, both studies using exactly the same observational data, viz., the 
surface brightness profile and 469 stellar radial velocities. 

• Parametric: a simultaneous fit of these radial velocities and of the surface 
brightness profile to a multi-mass King-Michie dynamical model provides mean 
estimates of the total mass for ui Cen equal to Mtot = 5.1 =t 0.6 x 1O 6 A^0, with 
a corresponding mean mass-to-light ratio M/Ly = 4.1 (Meylan etal. 1995). 

• Non-parametric: the potential and mass distribution infered in this method 
provide a total mass for a; Cen equal to Mtot = 2.9 ± 0.4 x 1O 6 A40, with a 
corresponding mean mass-to-light ratio M/Ly = 2.3 (Merritt etal. 1997). 

3. Universal Mass-to-Light Ratio? 

With a similar parametric approach applied to the observational constraints 
obtained for NGC 1835, an old Large Magellanic Cloud globular cluster, King- 
Michie models give Mtot = 1-0 ± 0.3 x 10 6 .M©, corresponding to a mean mass- 
to-light ratio M/Ly = 3.4 ± 1.0 (Meylan 1988, Dubath & Meylan 1990). 

These studies show that when the same kind of dynamical models (King- 
Michie) constrained by the same kind of observations (surface brightness profile 
and central value of the projected velocity dispersion) are applied to an old and 
bright Magellanic globular cluster, viz., NGC 1835, the results seem similar to 
those obtained in the case of Galactic globular clusters. Consequently, the rich 
old globular clusters in the Magellanic clouds could be quite similar in mass and 
M/Ly to the rich globular clusters in the Galaxy. 
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4. The Final Word from Thousands of Stellar Space Motions 

Recently, the amount of data related to stellar motions, viz., radial velocities 
and proper motions, has increased significantly. In a pioneering ground-based 
study of ijj Centauri, van Leeuwen et al. (2000) measured the individual proper 
motions of 7853 probable member stars, from photographic plates with epochs 
ranging from 1931 through 1935 and 1978 through 1983. An internal proper 
motion dispersion of 1.0 to 1.2 masyr -1 , equivalent to 25 to 29 kms -1 for a 
distance of 5.1 kpc, is found for members near the cluster center. This dispersion 
decreases to 0.3 masyr -1 , equivalent to 7.5 kms -1 in the outer regions. The full 
dynamical interpretation of these proper motions, combined with a few thousand 
radial velocities, is in preparation by the same group. 

Another group, using HST/WFPC2 images, has obtained slightly better 
measurements of proper motions for about 15,000 stars in the core of 47Tucanae, 
this within a time baseline of only 4 years, with 3 epochs separated by 2 years. 
See Anderson, King, & Meylan 1998 for a progress report after the second epoch. 
These data, combined with the radial velocities of about 5,000 stars, will pro- 
vide, as in the case of to Centauri, an insight into the dynamics of the core of 
47 Tucanae, with fundamental by-products such as cluster distance and photom- 
etry of variable stars and binaries. 

With their new proper-motion techniques and software (Anderson & King 
2000) applied to the HST/WFPC2 archive data for the first epoch and to their 
own observations for the second epoch, the members of the same group will soon 
have available similar sets of proper motions for about ten among of the richest 
and most nearby Galactic globular clusters. 

The above studies, with thousands of proper motions and radial velocities 
constraining dynamical models with three integrals of the motion as well as non- 
parametric ones, will allow a significant step forward in our understanding of 
the internal dynamics of massive star clusters. 



5. Mayall II = Gl, a Giant Globular Cluster in M31 

Mayall II = Gl is one of the brightest globular clusters belonging to M31, the 
Andromeda galaxy. Observations with HST/WFPC2 provide photometric data 
for the I vs. V — I and V vs. V — I color-magnitude diagrams. They reach stars 
with magnitudes fainter than V = 27 mag, with a well populated red horizontal 
branch at about V = 25.3 mag (Meylan etal. 2001). From model fitting, that 
study determines a rather high mean metallicity of [Fe/H] = - 0.95 ± 0.09, some- 
what similar to 47 Tucanae. In order to determine the true measurement errors, 
Meylan et al. (2001) have carried out artificial star experiments. They find a 
larger spread in V — I than can be explained by the measurement errors. They 
attribute this to an intrinsic metallicity dispersion among the stars of Gl, which 
may be the consequence of self-enrichment during the early stellar /dynamical 
evolutionary phases of this cluster. So far, only ui Centauri, the giant Galactic 
globular cluster, has been known to exhibit such an intrinsic metallicity disper- 
sion. This is, a phenomenon certainly related to the deep potential well of each 
of these two star clusters, which are massive enough to retain the gas expelled 
by the first generations of very massive stars. 
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The structural parameters of Gl are deduced from the same HST/WFPC2 
data. Its surface brightness profile provides its core radius r c = 0.14" = 0.52 
pc, its tidal radius r% — 54" = 200 pc, and its concentration c = log (rt/r c ) ~ 
2.5. Such a high concentration indicates the probable collapse of the core of 
Gl. KECK/HIRES observations provide the central velocity dispersion <t & s = 
25.1 kms -1 , with cr p (0) = 27.8 kms -1 once aperture corrected. 

Three estimates of the total mass of this globular cluster can be obtained. 
The King-model mass is Mr = 15 x 10 6 .M© with M/Ly ~ 7.5, and the Virial 
mass is Mvir = 7.3 x lO 6 ^© with M/L v ^ 3.6. The King-Michie model 
fitted simultaneously to the surface brightness profile and the central velocity 
dispersion value provides mass estimates ranging from Mkm = 14 x lO 6 ^© 
to 17 x W 6 M Q (Meylan etal. 2001). 

Table 1. Three mass determinations for Mayall II = Gl and wCen 



Mass Mayall II u> Centauri 



[W 6 M & ] [W 6 M Q 



King 15 4.3 

Virial 7.3 2.9 

King-Michie 13-18 5.1 



The spread between the three mass determination values listed in Table 1 
give a better idea than their individual formal (smaller) errors about their true 
uncertainties. The masses of both clusters are known to about a factor of two. 
Although not very precise, all of these mass estimates make Gl more than twice 
as massive as to Centauri, the most massive Galactic globular cluster. Gl is 
unique in M31 because of its projected location 40 kpc away from the center 
of the galaxy, but there are at least three other bright globular clusters in this 
galaxy which have velocity dispersions a j, s larger than 20 kms -1 , implying 
rather large masses. 

6. On the Origin of the Most Massive Globular Clusters 

Such large masses are related to the metallicity spread the origin of which is 
still unknown. It may come either (i) from metallicity self-enrichment in a mas- 
sive globular cluster, (ii) from primordial metallicity inhomogeneity in a binary 
proto-cluster cloud, followed by early merger, or (iii) from the fact that the 
present globular cluster is merely the remaining core of a previously larger en- 
tity, e.g., originally a dwarf galaxy subsequently pruned by dynamical evolution. 
Although uj Centauri is the best studied globular cluster, because of its size 
and relative proximity (~ 5.1 kpc), many conundrums remain: viz., (i) the 
metallicity spread among stars (Freeman & Norris 1981), (ii) its double Main 
Sequence (Anderson 1997), (iii) the different kinematics between metal-rich and 
-poor stars (Norris etal. 1997), and (iv) a correlation between metallicity and 
age, implying that this cluster enriched itself over a timescale ~ 3 Gyr (Hughes 
& Wallerstein 2000 and Hilker & Richtler 2000). 
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Table 2. Parameters for the massive globular cluster Mayall II = Gl 
and the nucleus of the dwarf elliptical galaxy NGC 205 (see Meylan 
etal. 2001) 



Parameters 


Mayall II 


NGC 205 


<7 P (0) 


27.8 kms" 1 


30 kms _i 


M v 


-10.94 mag 


-9.6 mag 


r c 


0.52 pc 


0.35 pc 


H(0,V) 


13.47 mag/arcs 2 


12.84 mag/arcs 2 
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Figure 1. Reproduction of Fig. 3 from Kormendy (1985) in which 
we have added the positions, using the four parameters in Table 2, of 
Mayall II = Gl (star) and of the nucleus of NGC 205 (triangle). These 
fall in all four panels right on the sequences defined by globular clusters 
and always away from the sequences defined by ellipticals, bulges, and 
dwarf ellipticals. 
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Let us consider (see Table 2 and Meylan et al. 2001) the following four pa- 
rameters relative to Gl: the central velocity dispersion <r p (0) = 28 kms -1 , the 
integrated absolute visual magnitude My = - 10.94 mag, the core radius r c = 
0.52 pc, and the central surface brightness /j,(0,V) = 13.47 mag arcsec -2 . The 
positions of Gl in the different diagrams defined by Kormendy (1985), using the 
above four parameters, always put it on the sequence defined by globular clus- 
ters, and definitely away from the other sequences defined by elliptical galaxies, 
bulges, and dwarf spheroidal galaxies (Fig. 1). The same is true for w Centauri 
(Meylan etal. 2001). 

Little is known about the positions, in these diagrams, of the nuclei of 
nucleated dwarf elliptical galaxies, which could be the progenitors of the most 
massive, if not all, globular clusters (Zinnecker 1988, Freeman 1993). The above 
four parameters are known only for the nucleus of one dwarf elliptical, viz., 
NGC 205, and their values put this object, in Kormendy's diagram, close to Gl, 
right on the sequence of globular clusters (see Table 2 and Fig. 1). This result 
does not prove by itself that all massive globular clusters are the remnant cores 
of nucleated dwarf galaxies. 

At the moment, only the anti-correlation of metallicity with age recently 
observed in to Centauri suggests that this cluster enriched itself over a time scale 
of about 3 Gyr (Hughes & Wallerstein 2000 and Hilker & Richtler 2000). This 
contradicts the general idea that all the stars in a globular cluster are coeval, and 
may favor the origin of u> Centauri as being the remaining core of a larger entity, 
e.g., of a former nucleated dwarf elliptical galaxy. In any case, by the mere fact 
that their large masses imply complicated stellar and dynamical evolutions, the 
very massive globular clusters may blur the former clear (or simplistic) difference 
between globular clusters and dwarf galaxies. 



7. Mass Estimates of Young Star Clusters 

So far, we have discussed only mass estimates related to old, rich star clusters. 
HST has triggered numerous studies of young, bright starburst clusters which 
may be quite massive. See, e.g., Holtzmann et al. (1992) in the case of NGC 1275 
clusters, and Schweizer & Seitzer (1993) and Whitmore etal. (1993) in the case 
NGC 7252. Conti & Vacca (1994) use HST/FOC UV imaging of various bright 
knots in the nuclear starburst region of the Wolf-Rayet galaxy H2 2-10. The 
luminosities of the knots are compared to those predicted from stellar population 
synthesis models for ages between 1 and 10 Myr, from which they estimate the 
masses of the knots to be between 10 5 and 10 6 Mq. In a more direct way, Boker 
etal. (1999) measure a velocity dispersion a b s = 33 ± 3 kms -1 in the nuclear 
star cluster of the face-on giant Scd spiral galaxy IC 342 and deduce, by fitting 
the central surface brightness profile, a cluster mass Mtot = 6 x 10 6 A4 Q . They 
infer a best-fitting cluster age in the range 10-60 x 10 6 yr. 

It is worth emphasizing that mass estimates of young star clusters are more 
difficult to obtain and more uncertain than those of older clusters. This is due 
to the frequent lack of velocity dispersion measurement and because of the large 
uncertainties inherent in the use of stellar population synthesis models for young 
stellar systems whose luminosities evolve quickly (see, e.g., Bruzual 2001). 
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8. The Slow Destruction of Star Clusters 

Numerical works have demonstrated that continual two-body relaxation within 
globular clusters combined with weak tidal encounters between globular clusters 
and the Galactic disk and/or bulge will lead, around each globular cluster, to 
the development of both a halo of unbound stars and tidal tails. All clusters 
observed, which do not suffer from strong observational biases, present tidal 
tails (Grillmair etal. 1995, Leon etal. 2000). These tidal tails exhibit projected 
directions preferentially aligned with the cluster orbit or towards the Galactic 
center, betraying their recent dynamical evolution through disk and/or bulge 
shocking. See also Odenkirchen etal. (2001). 

Recent theoretical work corroborate these observations. In N-body simu- 
lations of globular clusters moving in the Galactic potential well, Combes et al. 
(1999) observe that once the particles (stars) are unbound, they slowly drift 
along the globular cluster path and form two huge tidal tails. A cluster is al- 
ways surrounded by two giant tidal tails and debris, in permanence along its 
orbit. The length of these tidal tails is of the order 5 tidal radii or more. The 
orientation of these tidal tails is the signature of the last disk crossing and can 
constrain strongly the cluster orbit and the Galactic model. Each disk/bulge 
crossing may extract up to about 1% of the total mass of the star cluster, leading 
slowly but surely to its complete evaporation. The lighter the cluster and the 
stronger the tidal shocks, the faster the destruction process. 

9. Do Globular Clusters Have Dark Halo? 

Are globular clusters the most massive stellar systems without non-baryonic 
dark matter? Since both the velocity dispersion profile and the rotation curve 
of uj Centauri decrease towards the outer parts of this cluster, we may conclude 
that there is no dynamical evidence of any massive halo made of non-baryonic 
dark matter. 

In order to investigate further this point among the most poorly studied 
Galactic clusters, Cote et al. (2001) have obtained KECK/HIRES high-accuracy 
radial velocities of about 20 stars in six outer-halo globular clusters, located 
between 20 and 100 kpc from the center of our Galaxy. The velocity dispersions 
range between a b s = 1 and 5 kms -1 , with corresponding Ai/Ly values between 
1 and 4. This is exactly what is expected from the nearby well studied globular 
clusters such as to Centauri and 47 Tucanae. 

However there is one single, although conspicuous, exception in the sample 
of Cote etal. (2001): Pal 13 exhibits a velocity dispersion a b s = 2.6 ± 0.3 
kms -1 . With its low and uncertain total luminosity, Pal 13 has a corresponding 
Ai/Ly in the range 10 < A4/Ly < 40. This is quite unique for a globular cluster. 
Simulations show that such a result would be mimicked, with difficulties, only 
by an uncomfortably large fraction of binary stars. Binaries exist in globular 
clusters, but not in large quantities. We should apply Ockham's razor before 
invoking the presence of a massive dark halo. Such a high Ai/Ly value could be 
an indication of a cluster in the late phases of its dynamical evolution, when a 
large fraction of its total mass is made of white dwarfs, as predicted by Vesperini 
& Heggie (1997). It could also be explained by a velocity dispersion inflated 
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Figure 2. From bottom to top, a sequence of increasingly disturbed 
surface brightness profiles (two globular clusters and five dSph galaxies) 
corresponding to increasing M./Ly values. 



because of lack of Virial equilibrium: Pal 13 seems to be in the advanced stages 
of tidal disruption, a status supported by the recent observations of extra-tidal 
extension around this cluster (Siegel et al. 2001 and Cote et al. 2001). These two 
explanations do not exclude each other. 



10. Is Tidal Disruption an Ubiquitous Phenomenon? 

Leon et al. (2000) show that any Galactic globular cluster observed which does 
not suffer from strong observational biases (e.g., intervening absorption by Galac- 
tic cirrus along the line-of-sight) displays a pair of tidal tails. The tidal distur- 
bance may also be unveiled through its effect on the surface brightness profile. 
An isolated cluster suffering no tidal shock will have a normal King-model pro- 
file characterised by a core surrounded by an envelope with a steep profile. Any 
tidally-perturbed cluster will have an envelope profile departing from the King- 
model profile as if there were a very strong background of stars: the stronger 
the tidal effect, the higher the level of the background (i.e., the higher the de- 
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parture from King- model profile), which is made of stars escaping the cluster 
gravitational attraction. 

When observing a tidally perturbed globular cluster, there is a high prob- 
ability of measuring radial velocities of escaping stars, which are no more in 
Virial equilibrium with the cluster. Such stars would inflate the measured ve- 
locity dispersion. The stronger the tidal disturbance, the higher the inflated 
velocity dispersion, and the larger the mass-to-light ratio. This is exactly what 
is observed in Fig. 2, from Majewksi etal. (2001), which displays the surface 
brightness profiles of two globular clusters (NGC 288 with M/Ly = 2.1 and 
Pal 13 M/Ly ~ 20) and five dSph galaxies (Sculptor with M/Ly = 3.0, Leo I 
with M/Ly = 5.6, Leo II with M/Ly = 17, Carina with M/Ly = 31, and 
Ursa Major with M/Ly = 79). From the bottom panel (NGC 288) to the top 
one (Ursa Major) we observe a sequence of increasingly disturbed profiles, with 
higher and higher departures from the King models, with the most disturbed 
profile being for Ursa Major. Interestingly, this sequence of increasingly dis- 
turbed profiles correspond to a sequence of increasing mass-to-light ratios, as if 
the M/Ly values were a direct measure of the intensity of the tidal shocks. Is 
the evidence for dark matter in dSph galaxies evaporating with their stars? 
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Discussion 

Armandroff: You raised the possibility of tidal disruption leading to the large 
velocity dispersion in Ursa Minor. Models for the disruption of dSphs (Piatek & 
Pryor, 1995, AJ, 109, 1071; Oh etal., 1995, ApJ, 442, 142) predict a kinematic 
signature that resembles a rotation curve, as opposed to simply an increased 
dispersion. Large samples of radial velocities in Draco and Ursa Minor have not 
shown evidence for the kinematic signature of disruption. Do you have evidence 
for Ursa Minor or other dSphs matching the modeling? 

Meylan: The two papers you mention present results of N-body simulations 
which reach conclusions different from those obtained from other, more recent 
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N-body simulations (Kroupa, 1997, New Astron, 2, 139; Klessen &; Kroupa, 
1998, ApJ, 498, 143). Kroupa's work shows that dSphs are obtained under the 
extreme hypothesis that dSph progenitors are not dominated by dark matter 
but are significantly shaped by tides through many periastron passages; Pi- 
atek & Pryor simulated only one passage, and Kroupa's results are completely 
consistent with their findings for one passage only. There are now in Draco 
photometric observations of stars beyond its measured tidal boundary (Piatek 
et al., 2001 AJ, 121, 841). It will be essential to obtain, for a few of these dSphs, 
samples of thousands of stellar radial velocities. At the moment, I simply find 
the correlation between tidal disturbance intensity and A4/L values rather in- 
triguing. This definitely calls for more studies. 

Da Costa: Comment: Integrated spectrum of nucleus of NGC 205 is A-type, so 
its stellar population is very different from globular clusters like Gl. Question: 
What is the core radius of the NGC 205 nucleus and does if fit a King model? 
(different from say M32 with its power-law cusp). 

Meylan: About the comment: Yes, the stellar populations in the core of NGC 205 
are younger than those in Gl. This means that the position, in the panels de- 
fined by Kormendy (see Fig. 1 above), of a dynamical system older than about 
1 Gyr does not depend strongly on its age. Answer to the question: The core 
radius of the nucleus NGC 205 is 100 mas = 0.35 pc (Jones etal., 1996, ApJ, 
466, 742), similar to a very dense globular cluster easily fitted by a King model. 
The star cluster in the center of NGC 205 is very different from the nucleus of 
M32, which has, e.g., a much larger velocity dispersion a b s = 150 km s _1 (Joseph 
etal., 2001, ApJ, 550, 668). See also Meylan etal. (2001). 

Grillmair: Two comments: 1) The onset of tidal effects (e.g., break from normal 
profile to power-law tidal tail, velocity dispersion increase) can depend signifi- 
cantly on the orbital phase of the cluster (strongest effect at apogalacticon) . 2) 
Ben Moore (1996, ApJ, 461, L13) put an upper limit of M/L of NGC 7089 = 
M2, based on existence of tidal tails. 

Wenderoth: Around 15 years ago, the idea that uj Centauri is a merger was 
proposed. Since then, is there more evidence to support this idea or to reject it? 

Meylan: There was a paper by Icke & Alcaino, 1988, A&A, 204, 115, suggesting 
that uj Centauri could be the result of a merger, this in order to explain (i) its 
spread in metallicity, (ii) its strong flattening and (iii) its large mass. There 
is now more evidence about the intrinsic complexity of this globular cluster, as 
mentioned above in Section 6, although none of these observational facts sup- 
ports the merger scenario over the other two alternatives. If merger there was, 
it must have been very early in the life of the two merging proto-cluster clouds, 
since the mere encounter of two current globular clusters would not induce a 
merger, except in the case of orbital-parameter adjustment with vanishingly 
small probabilities. 



